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ABSTRACT 

We present a detailed X-ray study of the 2003 XMM-Newton observation of the High 
Mass X-ray Binary XTE J0421+560/ CI Cam. The continuum of the X-ray spectrum is 
well described by a flat power law (F = 1.0 ±0.2) with a large intrinsic absorbing col- 
umn (Nh = (4.4 ± 0.5) X 10^^ cm^^). We have decomposed the broad iron line into 
3 separate components: Fel-Ka, Fel-K/3 and FeXXlV-XXVKa. It is unclear how both 
neutral and almost fully ionised iron can exist simultaneously, however we suggest 
this could be evidence that the compact object is embedded in the circumstellar ma- 
terial. This doesn't appear to be consistent with the X-ray flux and spectrum of the 
source, which has remained essentially unchanged since the initial outburst. The iron 
abundance implied by the ratio of the neutral Fe-Ka and Fe-K/3 is compatible with 
solar. We search for lags in the neutral Fe-Ka with respect to the continuum and find 
marginal evidence for a lag at ~10 ks. We interpret this as the light crossing time of 
the torus which would suggest that the neutral iron is located at a radius of 10 AU. 
This result depends on several assumptions including the distance to the system, the 
inclination, the mass of system and the orbital period none of which are known with 
any great certainty. Better constraints on these system parameters and further obser- 
vations of this system are required to confirm this result. We discuss the nature of 
this system in light of our results and place it in context with other binary B[e] stars. 

Key words: X-rays: binaries, stars: emission-line. Be, techniques: miscellaneous, 
stars: individual: XTE J0421+560 



1 INTRODUCTION 

The X-ray transient XTE J0421+560 was discovered by the 
All-Sky M onitor on-bq ard RXTE during a strong outburst 
in 1998 by ISmith et"an It brightened to a peak flux of ^2 
Crab within a few hours before rapidly decaying, reaching 
quiescence in under two weeks. Observations at different 
wavelengths immediately after the outburst revealed simi- 
lar behaviour in the optical and radio bands with a peak 
flux density of almost 1 Jy at 8.30 GHz (e.g. IClark et al] 
12000!) . Initial radio observations suggeste d that there m ight 
be a superluminal expansion ( Hiellmin g et alj Il998h , but 
subsequent analysis showed that the radio emitting region 
was expanding with a veloci ty of only ~ 12,000 km s" 



( iMioduszewski & RupenI 20Qi) . Th e X-ray spectra se en by 
SAX jProntera et all Il998l) . RXTE fBelloni et al."l999^ and, 
more recently, XMM-Newton l lBoirin et aL .2002) are domi- 
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nated by an emission feature at 6-7 keV attributed to the Fe- 
Ka fluorescence line. 

The optical coimterpart of XTE J0421+560 is CI 
Camelopardali s (CI Cam ), a BO -2 supergiant B[e] star as 
classified by .Clark etalj ( |l999l) , iHynes et alj (2002) and 
iRobinson et alj ( I2OO2I) ," establishing this system as a High 
Mass X-ray Binary (HMXB). Common properties of stars ex- 
hibiting the B[e] phenomenon are the presence of forbidden 
emission lines in their optical spectra and a strong infra- 
red (IR) excess. T his IR excess is attrib uted to a hot circum- 
stellar dust shell jZickgraf etalj Il986l) . Until the launch of 
}NTEGRAL CI Cam was the only known sgB[e] compan- 
ion star in a HMXB. INTEGRAL has since i dentified another 
sgB[e]/XRB candidate, ICR J16318-4848 dCourvoisier et alj 
12003.) . 

The circumstellar gaseous environment of CI Cam 
has been imder mu ch scrutiny since the 1998 outburst. 
iRobinson et alj ( 12002) propose a two-component wind from 
the sgB[e] star to explain their observations; a cool, low veloc- 
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Table 1. XMM-Neioton EPIC observations of XTE J0421+560 on 2003 
February 24 



Camera 


Filter 


Read out 


Observation 


Exp. 






Mode 


Start End(UT) 


(ks) 


MOSl/2 


Medium 


FF, 2.6 s 


12:34 05:50 


59.5 


pn 


Thinl 


FF, 73 ms 


12:56 05:50 


60.3 



ity "iron" wind and a hot, high velocity wind. The low veloc- 
ity wind is dense, roughly spherical and extends to a radius 
of 13 to 50 AU. They note it is unclear as to how these winds 
with very different densities and velociti es co-exist. FoUow - 
ing a series of spectroscopic observations lHynes et all i2002l) 
also conclude that the emission lines must origin ate from 
several physically separate regions. Unlike Robins on et alj, 
they associate the H Balmer lines, He I and Fe II lines with 
the same region; a disk, or equatorial outflow region, viewed 
nearly pole on. This is consistent with the two component 
stellar wind model for supergiant B[e] stars proposed by 
IZickgrafetaL. (1985, 1986). 

Analysis of hi g h res olution optical spectra by 
iMiroshnichenko et alj l l2002h , observed almost 3 years 
after the outburst, revealed clearly resolved emission line 
profiles, displaying a triple peak structure. They suggest 
that this implies an intermediate inclination angle of the 
circumstellar disk with respect to the line of sight. Recent 
interferometric work in the n ear infrared ov e r an e ight year 
period from the outburst by iThureau et alj | |2009|) resolves 
the emitting region. They model this as a skewed elliptical 
Gaussian ring with semi-major axis a = 19.0 ± 1.2 AU 
(assuming a distance of 5 kpc) w ith an inclination angle 
of i ~ 67° . 'Barsukova et al.' teOOe") report a 19.41±0.02 day 
orbit period in the optical light curve of CI Cam, along with 
Doppler shifts in the Hell 4686 A line on the same time scale. 
The proposed period would put the semi-major axis of the 
binary orbit between 100-130 Rq: close to the surface of the 
star (based on mass and radii measurements of sgB[e] stars 
byrZickgraf et al. 1986). 

Unfortunately, distance estimates to XTE J0421+560 are 
uncertain with a wide range of values being adopted in the 
literature. This is primarily due to the lack of reliable spectral 
diagnostics available to determine the the luminosity class 
of the mass donor. Consequently, estimates range from 1- 
17 kpc. This in turn means that the X-ray luminosity of the 
outburst is not well constrained and as a result the nature 
of the compact object of this binary system is still unclear. 
We critically discuss the distance and nature of the system in 
Section 5. 

This paper reports on a 60 ks XMM-Newton observation 
of XTE J0421+560 during satellite revolution #0588 on 2003 
February 24. We discus s the dat a collected with th e EPIC 
MOS (Turner et ailUOOl ) and pn jStriider et 31112001 ') detec- 
tors and attempt to constrain the system geometry via spec- 
tral and temporal analysis. We also discuss the nature of the 
compact object and the system in light of these results and 
place it in context with other binary B[e] stars. 

2 OBSERVATIONS AND DATA REDUCTION 

Table [l] summarises the details of the EPIC observations. All 
detectors were operating in Full Window modes, giving a 



field of view of ~30' diameter. The data were processed using 
the XMM-Newton Science Analysis System v9.0 (SAS) along 
with software packages from FTOOLS v6.8. 

The MOS and pn observational data files were pro- 
cessed with emproc and epproc respectively. "Single and 
double" (PATTERN^4) pixel event patterns were selected 
for the pn detector. For the MOS instruments, "single" 
to "quadruple "(PATTERN^12) pixel events were selected. 
Source counts were extracted from a 30" radius region cen- 
tered on XTE J0421+560 and were compared with those 
from statistically identical background regions with extrac- 
tion radii 60" . Whilst several periods of high background ac- 
tivity are present in the data, the source is brighter than the 
background for all but the final ~13 ks of the observation. 
Hence we have not included this final section in the Ught 
curves reported here. The counts from the closest of the sta- 
tistically identical background regions (located on either a 
neighbouring CCD (pn) or the same CCD (MOSl and 2) as 
XTE J0421+560) were used as the final background regions 
for the light curves. The background subtraction was per- 
formed using the epiclccorr task which also corrects for bad 
pixels, vignetting and quantum efficiency. 

For the spectra, periods of high background activity 
were screened by removing any times when the count rate 
above 10 keV was >0.8 cts s^^ for the MOS detectors and 
>1.5 cts s^^ for the pn detector. Spectra were extracted us- 
ing a 24" extraction radius as recommended by the SAS task 
eregionanalyse, which maximises the signal to noise ratio 
based on the source counts. "Single and double" pixel events 
(PATTERN^4) were accepted for the pn detector, all bad pix- 
els and columns were disregarded (FLAG=0). For the MOS, 
"single" to "quadruple "(PATTERN^12) pixel events were 
selected with quality flag #XMMEA_EM. 

Background spectra were taken from regions with a 
60"radius away from the source on the same CCD. Sev- 
eral background regions were extracted and examined to en- 
sure they were statistically identical. The region from which 
the background spectrum for the pn detector was selected 
has a similar distance to the readout node as the source re- 
gion, as recommended bythe "EPIC status of calibration and 
data analysis" documen{j. For the MOS detectors, the same 
source region was used for both the MOSl and MOS2 de- 
tectors. The area of source and background regions were cal- 
culated using the backscal task. Response matrix files were 
created for each source using the task rmfgen and arfgen. 



3 SPECTRAL ANALYSIS 

The spectral analysis discussed in this work was performed 
using XSPEC ( Arnaud 1998) version 12.7.0. The spectra from 
the three EPIC detectors were fit simultaneously with the 
model parameters constrained to be identical across the three 
instruments. All models have an additional constant factor to 
account for the instrumental differences. All elemental abun- 
dances are set to the values of Wilms et al. (2000) unless oth- 
erwise stated. 

Figure [T] shows the spectrum of XTE J0421+560 both be- 
fore and after background subtraction along with the back- 

^ http: / /xmm2.esac.esa.int/ docs/ documents/ CAL-TN-0018.pdf 
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Figure 1. The 0.2-12.0 keV spectrum of XTE J0421+560. Top panel 
shows the background subtracted spectrum, bottom panel shows 
the source spectrum before background subtraction (in black) along 
with the spectrum of the backgroimd (red). All spectra are nor- 
malised. 

ground spectrum used. The spectrum of XTE J0421+560 is 
characterised by heavy absorption and a large emission fea- 
ture at ~6.5 keV.'Boirin et all {lOOf i report a soft (<2 keV) ex- 
cess in their 2001 observation. Whilst there does appear to be 
evidence for a soft excess in our background subtracted spec- 
trum, when the uncorrected source spectrum is compared 
with the backgroimd spectrimi, the shape of the spectrum 
and value of the data points below 3 keV are almost identi- 
cal. A visual inspection of the X-ray image shows no detec- 
tion of the source below 3 keV. Whilst we cannot rule out the 
existence of a soft excess in the spectrum of XTE J0421+560, 
we believe that this could be an artifact of the backgroimd 
subtraction and so do not include the <3 keV data points in 
our model fits. 

The continuum of the spectrum was first fit by exclud- 
ing the data around the iron line region. The absorption was 
modeled with two elements, a Galactic foreground compo- 
nent, NH,Gai fixed to 4.5 x 10^^ cm~^ jPickey & Lockmaiil 
Il990h and a separate column density, NH,i, intrinsic to the 
source which was allowed to vary. A simple absorbed power- 
law model {phabs*phabs*powerlaw in XSPEC) yielded a re- 
spectable fit with a value of 135.6 for 128 degrees of free- 
dom (dof; Xr = 1-06). A partially covered power law model 
was also fit to the data {phabs*pcfabs*poiverlaw in XSPEC). In 
this model the underlying powerlaw is fractionally covered 
by an absorbing material with the uncovered emission ab- 
sorbed by the interstellar medium. This did not improve the 
continuum fit with a value of 144.4 for 127 dof (Xr = 1-14). 
We also fit the continuum with a power law partially cov- 
ered by an ionised absorber (phabs*pcfabs*zxipcf in XSPEC). 
This produces a statistically better fit than the partial cover- 
ing model alone with a x^ of 138.4 for 133 dof (xr = 1-04), 
but with unfeasible model parameters (F — 3.4l} g). Again, 
the covering fraction and ionisation parameter both suggest 
that an ionised partial covering model is imnecessary. 

The region around the iron line was then included in the 
model fit. A single (unabsorbed), intrinsically narrow, Gaus- 
sian was added to the model with the energy of the line 
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Figure 2. The EPIC spectrum of XTE 10421+560 around the iron line 
along with the single and triple Gaussian models fit to the data. 



allowed to vary. The resulting fit has a Ka line at 6.43 ± 
0.02 keV. The fit has a x^ of 196.7 (xr of 1.17 for 168 dof). 
A better fit is produced by allowing the Gaussian width (x 
of 167.3; Xr of 1.00 for 167 dof), with a broad iron line at 
6.48 ± 0.04 with a = 0.23^0.05. Figure |2] shows the iron line 
region along with the best fit to the model and the broad iron 
line. Whilst statistically better than a single, narrow Gaus- 
sian, this fit has significant residuals and appears to be fit- 
ting the iron edge at 7.1 keV within the Gaussian. We do not 
believe this fit is physical. 

A second intrinsically narrow Gaussian was added to 
the model at 6.67 keV, consistent with ionised Fe-Ka. This 
further reduced the x' to 171.0 (xr=l-03 for 166 dof). Finally, 
a K/3 line was added to the model fixed at 7.06 keV. The 
normalisation was initially constrained to be 0.13 of that of 
the Ka line (Kaastra & Mewe 1993). This reduced the x^ to 
164.7 (xr=0.98 for 168 dof). When the normalisation and en- 
ergy are imconstrained the line energy remains essentially 
unchanged {EKfj = T.OSIq q^) and the ratio between K/? and 
Ka rises to 0.3±0.2, consistent with the theoretical value. 

Determining the signifi cance of the iron l ines in the spec- 
trum is not straightforward. [Protassov et al.l l l2002h in partic- 
ular warn against using the F-test to determine the signifi- 
cance of an emission line. Allowing the line energy to be free 
(as in our fitting) means that the models are no-longer nested 
and so the energy of the line cannot be incorporated as a free- 
parameter in the F-test. Instead, we need to allow for the fact 
that we have conducted a number of trials, A'^, over the en- 
ergy range of the line. Conservatively, we estimate that the 
chance detection p robability is at most ~ 1 — (1 — P)^ (see 
iPorquet et al J 12004 for a full discussion). Since the line en- 
ergies are within the instrument resolution of the expected 
values of 6.4 keV and 6.7 keV and not simply at some arbi- 
trary redshift, the number of independent trials that we have 
effectively searched over is small. Here we adopt a value of 
3 for iV. 

The data were fit with the simple absorbed power law 
model and the narrow iron lines added incrementally, as de- 
scribed above. They were then frozen at the best fit energy 
from the phabs*((phabs*po)+Gaus+Gaus+Gaus) model and an 
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Table 2. Best fit parameters for the phabs*((phabs*po)+Gaus+Gaus+Gaus) model fits to the full 3.0-12.0 keV spectrum. Errors, where reported, are 
the 90% confidence level. 
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Flux 
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171.0/166 


(4.4 ±0.5) X 10^3 


1.0 ±0.2 


6 41+0-03 


(6 ± 1) X 10-0 


6 67+0-02 
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7 05+0.08 
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Table 3. Estimates of the significance of the iron lines based on the 
two methods described in the text. Errors, where reported, are the 
90% confidence level. 
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Figure 3. EPlC-pn (black), EPIC-MOSI (red) and EFIC-MOS2 (green) 
spectra, along with best t phabs*((phabs*po)+Gaus+Gatis+Gatis) model 
to the data (solid line). Broken line shows the absorbed power law 
component. 



F-test was performed after the addition of each line. As such 
each model is compared with the previous, simpler model. 
We also estimate the significance using the equivalent widths 
of the lines and their errors. The results of both of these meth- 
ods are listed in Table|3] Both the neutral and ionised Ka lines 
appear to be clearly significant. There appears to be only a 
marginal detection of the expected K/? line, nonetheless the 
flux of the line is consistent with Solar metallicity. 

The continuum and line parameters for the best fit 
model are listed in Table|2] Figure |3] shows the best fit model 
to the spectrum. The final model fit has a photon index of 
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Figure 4. Background subtracted 3.0-12.0 keV light curve of 
XTE J0421+560 from the EPlC-pn detector The broken line shows 
the median count rate of the light curve. 



r = 1.0±0.2 and an intrinsic absorption of Nh = 4.4 ± 0.5 x 
10^3 cm-2, with a of 161.6 (Xr=0.98 for 164 dof). 

We als o attempted to fit the thermal model of 
llshida et alj 11504), derived from an ASCA observation, to 
the data. In this model emission from two regions of hot 
diffuse gas are photoelectrically absorbed independently 
(phabs*mekal+ph abs*mekal in XSPEC) wi th the elemental 
abundances of Anders & Grevesse 11989'). The parameters 
were frozen to the best fit values of Ishida et al. ( 2004), with 
only the normalisation allowed to vary, and fit to the 0.5- 
10.0 keV spectrum. This allowed for a direct comparison. The 
resulting fit was poor with a of 1023.6 (xr=5.62 for 182 
dof), rising to 1040.6 ( x?=6.09 for 171 dof) when only the 
3.0-12.0 keV energy range was considered. Allowing the ab- 
sorption and plasma temperature to vary improved the fit 
(X^=308.2; Xr=l-7 for 178 dof) but with unfeasible model pa- 
rameters (e.g. Nh ~ lO^'' cm-^, kT~80 keV). 



4 TIMING ANALYSIS 

Figure H] shows the background subtracted 3.0-12.0 keV 
light curve for XTE J0421+560, with a bin time of 500 s, 
from the EPIC-pn detector. The total net coimts for the hght 
curve, after background subtraction and excluding triple and 
quadruple pixel events, is ~5000 coimts. This corresponds 
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Figure 5. light curves of the background subtracted low energy con- 
tinuum (4.5-6.0 keV), the Ve-Ka line (6.2-6.55 keV) and the high en- 
ergy continuum (7.5-9.0 keV). 



to an average net count rate of 0.107±0.002 counts s ^, sub- 
stantially brighter than the 2001 observat ion which found 
0.024±0.002 counts s"^ jBoirin et al.ll2002h . The Ught curve 
has a fractional root-mean-squared v ariability amplitude of 
0.37 (calculated using equation (10) of lVaughan et alj|2003h . 

Figure|5] shows the source light curve split into three dif- 
ferent energy bands. The top panel shows the light curve for 
the lower energy continuum (4.5-6.0 keV), the middle panel 
shows the light curve in the range we have defined as the 
Fe-Ka line from our spectral fits (6.20-6.55 keV) and the bot- 
tom panel shows the light curve for the higher energy con- 
tinuum (7.5-9.0 keV). The Fe-Ko line spans a much narrower 
range of energies than the continuum light curves. This is to 
avoid contamination from the hypothesised ionised Ka line 
and the K/3 line, as well as to minimise any contribution of 
the underlying continuum to the cross correlation function 
(CCF). Above ~9.0 keV the signal to noise of the system re- 
duces and so these data are not included in our analysis. 

Cross-correlation analysis was performed on the three 
light curves mentioned above and an additional low energy 
continuum light curve ranging from 3.0^.5 keV. The IDL 
routine c_correlate.pro from the IDL Astronomy User's Li- 
brar}3 was used throughout. The continuum light curves 
were correlated with the sum of the other 3 continuum light 
curves (i.e. the rest of the spectrum without the iron line com- 
plex). The Fe-Ko line light curve was correlated with the sum 
of all the continuum light curves. Figure|6]shows the CCFs of 
the light curves with respect to the continua. The solid black 
lines are the same CCFs, smoothed using a moving average 
method in an attempt to minimise the noise. A positive time 
delay corresponds to a lag in the variation of the light curve 
with reference to the continua. The Fe-Ka CCF appears to be 
more asymmetric than its continuum counterparts. More ex- 
plicitly, the minimum present in the continua CCFs around 
^10 ks does not seem to be present. 

Whilst the energy range of the Fe-Ka line has been delib- 
erately defined to keep continuum contamination at a mini- 
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mum, inevitably there will be some contamination as there 
is no way to separate the line emission from the continuum 
emission at these energies. As such we expect to see some 
correlation of the light curves at zero lag. However, an asym- 
metr y in the CCF is ofte n interpreted as evidence for a lag 
(e.g. Pete rson et alj |l998) and so we explore the possibility 
that the shape of the CCF could be due a transfer function 
which includes some power at positive lags. 

We define the "asymmetry factor", Asym, of the CCFs to 
be the ratio of the difference of the half widths of the peak at 
tenth maximum, to the sum of the half widths. This gives a 
number between -1 and 1, with zero being perfectly symmet- 
rical. A positive value of Asym indicates that the gradient of 
the peak is steeper on the negative lag side and vice versa for 
a negative value. The Asym values of the continua CCFs are 
all very similar, ranging from -0.04 to -0.08, suggesting they 
are essentially symmetrical. However, the Asym value for the 
Fe-Ka CCF is 0.39, i ndicating a m oderate level of asymmetry. 

The method of I ZhangI (12002') was then used to see if the 
asymmetry seen in our Fe-Ka CCF could be caused by Pois- 
son noise. Monte-Carlo simulations were performed to as- 
sess the significance of this asymmetry. Ten thousand red- 
noise light curves were ge nerated with a po wer law slope 
of -2.0 using the method of 'Timmer & Koenie 1 1995) and the 
IDL routine rndpwrlc.pro. The light curves were initially sim- 
ulated with a duration ten times longer than that of the ac- 
tual data and were then cut down to the observed duration 
to account for red-noise leakage. The simulated light curve 
was then used to create four light curves with the statisti- 
cal properties of the observed light curves, resulting in four 
identical light curves distinguishable only by their mean and 
standard deviation. Noise was then added to each data point 
of the simulated light curves by adding a Gaussian devia- 
tion with a mean of zero and a variance equal to that of the 
corresponding data point in the real light curve (similar to 
the method of lUttley et al.ll2003h . Any negative data points 
were then set to zero. The three continuum light curves were 
then combined and cross-correlated with the fake Fe-Ka line 
light curve in the exact same manner as the observed data. 
The asymmetry factor and peak value of the smoothed CCFs 
were then calculated. 

Of the lO'* CCFs generated in the Monte Carlo simula- 
tion, 303 were observed to have Asym > 0.39 suggesting a 
significance of ~97% or 2.1(7. This result remains at a similar 
significance (2.3ct) when a red noise power law index of -1.0 
is considered. The asymmetry factors and peak values of the 
continuum CCFs all fall within the central 50% of the distri- 
bution. Whilst this result is not formally significant, the low 
signal to noise of our data makes a convincing result imlikely. 



5 DISCUSSION 

The continuum parameters of ou r best fit model ar e con- 
sistent with those reported by iBoirin et alj ( |2002|) . The 

^+0.04\ 



-0.05 -I 



total 3.0-10.0 keV flux of our model is (1.36 
10~^^ ergs cm~^ s~^ corresponding to an X-ray luminosity of 
Lx(3~iokeV) = (4.i1q;2) X 10^^ ergss"^ atadistanceof 5kpc. 
This is a slight increase on that seen in the previou s XMM ob- 
servation (Lx(o.5-iokeV) = 3.5 x 10^^ ergs s"^; Boir in et all 
I2OO2 ), however as there are no errors quoted on this value we 
cannot say whether this increase is significant. Whilst these 
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Figure 6. The cross-correlation functions of the light curves. The continua light curves (3.0-4.5 keV, 4.5-6.0 keV and 7.5-9.0 keV) were cross- 
correlated with the sum of the remaining continua light curves whilst the Fe-Ka (6.2-6.55 keV) light curve was cross-correlated with the sum of 
all the continua light curves. The soUd black line is smoothed CCF. 

luminosities do not cover the same energy range, as the sys- 
tem is undetected below 3 keV in our XMM observation it 
is imlikely that the I'x(o.5-iofceV) will vary much from the 
luminosity reported above. The flux values used to calcu- 
late the luminosities are derived from the spectral fits and so 



are also model dependant. The light curve seems to contra- 
dict the results from the spectral fits to the data, suggesting 
a factor of 4.46±0.09 increase in the average count rate over 
the two year period between data sets. A BeppoSAX obser- 
vation in 2000 put an upper limit on the X-ray luminosity 
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Viewed from Side: 




Viewed down Pole 




Figure 7. Sketch of the proposed geometry for CI Cam. The top figure 
shows the cross-section of the torus, whist the bottom figure shows 
the system as viewed down the pole. The thicker black ring at radius, 
r, is the assumed location of the neutral iron. The central object in this 
case is the sgB[e] star The position of the compact object is believed 
to be clos e to the surface of the o ptical star and thus lying well within 
the torus <Barsukova et alj2006l) . 



of L 



X{l-10keV) 



<2.5 X 10^=* ergs s"^ 



jParmar et ai]l200(]h . As 
such, the X-ray emission from the source does appear to have 
shown a moderate level of variability over the last few years. 

The improved signal to no ise has allowed the apparently 
broad Gaussian, reported by iBoirin et al.l b002l) , to be de- 
composed into three intrinsically narrow components, neu- 
tral Fe-Ka and K/3 and helium like FeXXV Ka. It is unclear 
how both neutral and almost fully ionised iron can co-exist 
and only seems possible if the ionised line emission arises 
from a physically separate region from the neutral emission, 
possibly local to the compact object. This would require the 
compact object to be in the same plane as the equatorial out- 
flows of the companion stars and seems inconsistent with the 
X-ray activity of the system. The iron abundance implied by 
the ratio of the neutral Ka and K/3 line fluxes is just consistent 
with solar within errors. 

Figure [7| shows t he current hyp othesised geometry for 
CI Cam proposed by iThureau et alj (2009), consisting of a 
binary system well within the torus. If we assume that the 
high column density originates from the leading (i.e. clos- 
est) edge of the surrounding torus (see Fig. then the in- 
clination, i ^ 67°, and semi-major axis, a=19 ± 1.2 AU 
iThureau et al.ll200^ . give a lower limit to the height of the 



torus of 8.1±0.5 AU. This is highly dependent on both dis- 
tance and inclination angle, neither of which are secure. 

The Kq fluorescence line of iron is already considered 
a powerful diagnostic tool for probing the central regions of 
both Gala ctic black hole candidate s ystems and active galac- 
tic nuclei ^Reynolds & Nowa5l2003l) . A technique known as 
reverberation mapping involves searching for time depen- 
dent changes in the continuum, emitted by the primary X-ray 
source, and then looking for the echo, or "reverberation", in 
the emission lines located further out. The relationship be- 
tween the continuum light curve and the emission line light 
curve is characterised by a geometry dependent "transfer 
function". Some idea of the spatial extent of the emission line 
region can be achieved by cross-correlating the light curves 
of the emission line, L{t + At), and the continuum, C{t). 
The function produced wiU be maximised at the lag, Ar, be- 
tween the emission line and the X-ray continuum. If we as- 
sume that all the neutral iron is located in the torus at radius 
r, that the compact object is in the centre of the torus (i.e. that 
the semi-major axis of the orbit is negligible compared to the 
semi-major axis of the torus) and that the X-ray emission is 
isotropic, then the range of lags seen is simply 



Ar = -(1 - cos6») 



(1) 



where 9 is the angle, between the line of sight and the X-rays 
incident on the torus (see Fig0. A more in depth descrip- 
tion along with the fundamental princip les and assumptions 
of reverberation mapping is given by iBlandford & McKei 
(1982) and Peterson (1993). 

There appears to be evidence for an asymmetry in the 
CCF of iron Kq energy range with the continuum, however 
much higher signal to noise data is needed to confirm this. 
If we consider the light curve of the Fe-Ka energy range, 
L{t + At), as the the sum of two lightcurves, Li {t + At) rep- 
resenting the line emission and L2{t) representing the con- 
tinuum "underneath" the iron line in the same energy range, 
then the CCF produced can be expressed as the sum of two 
CCFs: 



i: 

r 
/ 

J — ( 



L(i + At) C(t) dt 



{Li{t + At) + L2{t)) C{t) dt 



Li{t + AT)C{t)dt + 



f 



L2(t)C{t)dt (2) 



The first component of equation |2] represents the un- 
normalised CCF of the pure line emission (i.e. the CCF is not 
normalised by the root mean square of the light curve in each 
band) and the second component the un-normalised CCF of 
the underlying continuum. If this asymmetry is real and the 
continuum behaviour is simple and well represented by the 
CCF of the 7.5-9.0 keV light curve with the continuum (cho- 
sen as it is these photons which will be directly promoting 
the K-shell electrons into the L-shell), then subtracting the 
un-normalised 7.5-9.0 keV CCF from the the im-normalised 
Fe-Ka CCF will leave the CCF of the pure Fe-Ka emission 
with the continuum. 

Figure |8] shows the result of subtracting the un- 
normalised 7.5-9.0 keV CCF from the the un-normalised Fe- 
Ka CCF (consistent with the mathematical proof above). The 
poor signal to noise is in part a resul t of the 0.34 fluorescence 
yield of iron iBambynek et alj[l972l) : only ~1 in 3 of the in- 
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Figure 8. The CCF produced when the un-normalised 7.5-9.0 keV 
CCF is subtracted from the un-normalised Fe-Ko CCF. 

cident photons with energies >6.4 keV will result in a Ka 
photon. However, there does appear to be a distinct peak at 
5-10 ks. 

The most likely values of r from equation[T]are ~ and 
~ 2r, where r is the radius of the torus. If the peak is in- 
terpreted as the lag at 2r then this would indicate that the 
neutral iron is located at a radius of ~10 AU from the cen- 
tre of the torus. This result relies on several assumptions and 
parameters which are by no means secure. T he derived semi- 
majo r axis of the orbit (from the period of iBarsukova et all 
l2006 h depends on the cubed root of total mass of the system. 
As such large changes in the total mass of the system only 
produce small alterations in the semi-major axes. If we con- 
sider a companion star of mass 50 ± 30 Mq and a compact 
object with a mass 5 ± 5 Mq (consistent with observed neu- 
tron star and stellar mass black hole mass estimates), then the 
semi-major axis is 0.5 ±0.1 AU. This is a small fraction of the 
semi-major axis of the torus and so the approximation that 
the compact object is at the centre of torus seems valid. How- 
ever, in reality the neutral iron is unlikely to be located in 
an infinitesimally thin ring in the torus, it seems more likely 
that there exists an inner radius at which neutral iron can 
form. In this case a lag would also have a more complex den- 
sity dependence than the one that is considered here. Both 
these factors will serve to broaden the peak seen at 2r, and 
so the derived radius is likely the radius where the neutral 
iron density is at a maximum. 

5.1 The nature of CI Cam 

A significant obstacle to understanding the nature and be- 
haviour of CI Cam is the uncertaint y in the distance, with 
estimates r anging from ~ 11 k pc iBarsukova etal.ll2006h 
to >10 kpc (Robinson et alj|2002[) . The consequential uncer- 
tainty in the luminosity impacts on the nature of both the 
mass donor and the accretor, with the former authors favour- 
ing a B4III-V + white dwarf model and the latter a sgB[e] 
star + relativistic accretor (neutro n star or black hole ). In sup- 
port of the former interpretation, llshida et alj J2004') infer the 
presence of a white dwarf from analysis of the ASCA X-ray 
spectrum. This was fit with an optically thin, thermal hard X- 



ray model, characteristic of cataclysmic variables which con- 
tain a white dwarf. Our analysis (Section 3) clearly favours 
the power law model, with similar parameters to those 
reported for the other known sgBfej HMXB, ICR J16318- 
4848 (Matt & Guainazzi |2003|; iBarragan et"!!] |2009|) . which 
Ipilliatre & Chatyl 120041) claim hosts a neutron star based on 
the ratio of the X-ray and radio flux (though they note that 
the xmique environment means they are cautious about ap- 
plying such a relationship). The photon index reported is also 
consiste nt with thos e found in neutron star Be /X-ray bina- 
ries by lHaberl et alj (0.6-1.4; 2008). 

5.2 Comparison to other B[e] supergiants 

To date only two other supergiant B[e] stars have been de- 
tected as X-ray sources; the HMXB ICR J1 6318-4848 (e.g. 
iFiUiatre & Chatyl |2(]04'; 'Chatv & Rahoui"2012) and the clus- 
ter member Wdl-9Tskinner et al. 2006; Clark et al. 2008). ICR 
J16318-4848 shows a striking similarity to CI Cam, with a 
quiescent X-ray spectra best fitted with an absorbed, trun- 
cated power law and a strong Fe-Ka line. There is no evi- 
dence for the ionised iro n line in the higher resolution Suzaku 
data of IGR J1631-4848 I Barragan et al. 2009) suggesting that 
this could be unique to CI Cam. The X-ray spectrum of Wdl- 
9 instead demonstrates a spectrum consistent with emis- 
sion from an op tically thin thermal plasma (kT ~ 2.3keV; 
IClarketalJl200g ). This is suggestive of emission from a col- 
liding wind binary rather an accreting binary, consequently 
we do not discuss this source further. 

Uncertainties in the distances to both CI Cam and IGR 
J16318-4848 mean that it is difficult to compare quiescent lu- 
minosities, but our results suggest that they appear broadly 
comparable (to within a factor of a few) under the assump- 
tion that CI Cam is located at ~5 kpc and IGR J16318-4848 
is wit h in the 1.6-4 kpc range f avoured by IChaty & Rahouil 
i2012h iMatt & Guainazzil |2003|. Sect. 3) implying that the 
quiescent fluxes in both CI Cam and IGR J16318-4848 arise 
via a similar mechanism. Both sources have undergone X- 
ray "flares" with current temporal sampling suggesting that 
these events are rare and aperiodic; only two flares sepa- 
rated by ~ 9 y rs have been associated with IGR J16318-4848 
(1994 and 2003; ICourvoisier et aI.l2003l:lMurakami et al.l2003) 
and just the one event in 1998 with CI Cam (noting that 
photometric observations betwe en 1989-92 also s how no ev- 
idence for flaring in this period: IClark et aill2000h . Moreover 
the peak lumi nosity also differs between the systems; with 
iFilliat re & Chatv suggesting Lx ~ 2.6 x 10^" ergs s^^ for 
IGR J163818-4848 (assuming a distance o f 4 kpc) compared 
to Lx ~ 3 X 10^** ergs s~^ for CI Cam (e.g. lHynes et alj2002h . 

Both CI Cam and IGR J16318-4848 appear to be rather lu- 
minous supergiants, with comparable bolometric luminosi- 
ties. Both sources also support a mid-IR excess attributed 
to circumstellar gas. Although different modeling assump- 
tions (optically thin/ thick emission, spherical or disc geome- 
try and composition) make direct comparison of the physical 
properties (e.g. inner radius, dust mass) of the dusty com- 
ponents difficult. The dust temperature at the inner radius 
of the envelope appears rath er high in bo th cases (~ 767K 
for IGR T1 6318-4848 iChaty & Rahouill2012l) and ~ 1550K for 
CI Cam i Thureau et al.i i2009)). The geometry of the dusty 
en velope has bee n xmam biguously determined for CI Cam 
by iThureau et alj b009h , who showed that it resides in a 
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ring or torus with irmer radius ~ 8 AU (at 5 kpc). No con- 
clusions as to the dust geometry in IGR J16318-4848 were 
drawn by Moon et al. (2007), although by analogy to other 
B[e] stars iChaty & Rahoui. (i2012i) adopt a disc like geometry 
wit h inner radius ^ 1 AU /kpc. If the ~19.3 day orbital pe- 
riod Barsuko va et al. U2OO6') propose for CI Cam is confirmed, 
th e binary would also lie interior to the dusty disc mapped 
by IXhureau et alj |2009). Unfortunately no orbital period is 
known for IGR J16318-4848. 



6 CONCLUSIONS 

We have presented a detailed analysis of the spectral and 
timing properties of XTE J0421+560/CI Cam. The contin- 
uum of the spectrum is well described by a heavily absorbed 
{Nh = (4.4 ± 0.5) X lO^^cm"^ power law with photon in- 
dex consistent with those seen in neutron star Be /X-ray bi- 
naries (r = 1. ± 0.2). We also a ttempted to fit the white 
dwarf model of llshida et al.l j2004 l) but could not achieve an 
acceptable fit or sensible model parameters. The apparently 
broad iron line has been decomposed into three, intrinsically 
narrow lines, neutral Fe-Ka and K/3 and almost completely 
ionised FeXXV Ka, suggesting multiple emitting regions. The 
quiescent spectrum seems to point to a neutron star or black 
hole accretor, however this is inconsistent with the X-ray ac- 
tivity of the system and we cannot suggest a mechanism for 
the large X-ray flare seen in 1998. 

Further observations of XTE J0421+560/CI Cam of in- 
creased duration are required to confirm or refute the pres- 
ence of the lag suggested in this paper. Better quality data 
could be used to solve the transfer equation and could be 
used to explore the possibility of a lag between the neutral 
and ionised Fe-Ka emission, both of which are beyond the 
scope of this work. Better constraints on other system param- 
eters, such as distance, inclination and system mass, would 
also help to constrain the geometry of this unusual system 
using this method. 
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